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Abstract 

In the framework of Monitoring by Ocean Drifters (MONDO) Project, 
a set of Lagrangian drifters were released in proximity of the Brazil Cur- 
rent, the western branch of the Subtropical Gyre in the South Atlantic 
Ocean. The experimental strategy of deploying part of the buoys in clus- 
ters offers the opportunity to examine relative dispersion on a wide range 
of scales. Adopting a dynamical systems approach, we focus our atten- 
tion on scale-dependent indicators, like the finite-scale Lyapunov exponent 
(FSLE) and the finite-scale (mean square) relative velocity (FSRV) be- 
tween two drifters as function of their separation, and compare them with 
classic time-dependent statistical quantities like the mean square relative 
displacement between two drifters and the effective diffusivity as functions 
of the time lag from the release. We find that, dependently on the given 
observable, the quasigeostrophic turbulence scenario is overall compati- 
ble with our data analysis, with discrepancies from the expected behavior 
of 2D turbulent trajectories likely to be ascribed to the non stationary 
and non homogeneous characteristics of the flow, as well as to possible 
ageostrophic effects. Submesoscale features of ~ 0(1) km are considered 
to play a role, to some extent, in determining the properties of relative 
dispersion as well as the shape of the energy spectrum. We present, also, 
numerical simulations of an OGCM of the South Atlantic, and discuss 
the comparison between experimental and model data about mesoscale 
dispersion. 
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1 Introduction 



Detailed investigation of Geophysical flows involves experimental campaigns in 
which buoys, in the ocean, or balloons, in the atmosphere, are released in order 
to collect Lagrangian data against which theories and models can be tested. 
Questions concerning oil spill fate, fish larvae distribution or search and res- 
cue operations are only a few examples that make the study of advection and 
diffusion properties not only a challenging scientific task, but also a matter of 
general interest. 

In the past years, an amount of Lagrangian data about the South Atlantic 
Ocean (SAO) was collected thanks to the First Global Atmospheric Research 
Program (GARP) Global Experiment (FGGE) drifters, released following the 
major shipping lines, the Southern Ocean Studies (SOS) drifters, deployed in 
the Brazil-Malvinas Confluence (BMC) and the Programa Nacional de Boias 
(PNBOIA) drifters [Brazilian contribution to the Global Oceans Observing Sys- 
tem (GOOS)], released in the Southeastern Brazilian Bight (SBB). These data 
allowed estimates of Eddy Kinetic Energy (EKE), integral time scales and diffu- 
sivities (Piola et al. 1987; Figueroa and Olson 1989; Schafer and Krauss 1995). 
Despite the relatively uniform coverage, the boundary currents resulted poorly 
populated by buoys; furthermore, all previous studies about drifters in the South 
Atlantic have concerned one-particle statistics only. In this regard, in the frame- 
work of Monitoring by Ocean Drifters (MONDO) Project, a recent Lagrangian 
experiment, consisting in the release of a set of 39 World Ocean Circulation 
Experiment (WOCE) Surface Velocity Program (SVP) drifters, was planned in 
relationship with an oil drilling operation in proximity of the coast of Brazil, 
around (24°S, 44°W). Part of the drifters were deployed in 5-element clusters, 
some of them with initial drifter separations smaller than 1 km. This set of 
satellite-tracked Lagrangian trajectories offers, now, the opportunity to revisit 
advective and diffusive properties characterizing the current systems explored 
by the drifters. From the analysis of trajectory pair dispersion we can extract, 
in principle, information about the dominant physical mechanism acting at a 
certain scale of motion (e.g. chaotic advection, turbulence, diffusion). A thor- 
ough description of the oceanography of the South Atlantic Ocean, particularly 
of the main circulation patterns and of the mass transport properties, can be 
found in Peterson and Stramma (1991); Campos et al. (1995); Stramma and 
England (1999). The major feature characterizing the central region of the SAO 
is the large anticyclonic (anticlockwise) circulation known as Subtropical Gyre 
(SG). Other relevant surface current systems are: South Equatorial Current 
(SEC), Brazil Current (BC), Malvinas Current (MC), South Atlantic Current 
(SAC) and Benguela Current (BgC), as shown in Fig. [T] 

In the quasigeostrophic (QG) approximation, valid for relative vorticities 
much smaller than the ambient vorticity because of the earth's rotation, some 
theoretical arguments would predict that, from the scale of the forcing at which 
eddies are efficiently generated by instability, e.g. the Rossby radius of deforma- 
tion, both a down-scale enstrophy cascade and an up-scale energy cascade take 
place, corresponding to energy spectra E(k) ~ fc -3 and E(k) ~ fc -5 / 3 , respec- 
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Figure 1: Large-scale surface circulation in SAO (courtesy of Integrated Ocean 
Drilling Program; available online at |http: / /www.iodp.org/| . 
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Figure 2: Overall view of the trajectories of MONDO drifters. 

tively (Kraichnan 1967; Charney 1971). From a relative dispersion perspective, 
in the forward cascade range, the mean square relative displacement between 
two trajectories grows exponentially fast in time (Lin 1972), whereas, in the 
inverse cascade range, it follows the ~ t 3, power law (Obhukov 1941; Batchelor 
1950). Possible deviations from this ideal picture may reasonably come from 
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the non homogeneous and non stationary characteristics of the velocity field: for 
example, in the case of boundary currents, as well as from ageostrophic effects. 
At this regard, one presently debated issue is the role of submesoscale vortices 
(McWilliams 1985) [velocity field features of size ~ O(l) km] in determining 
the shape of the energy spectrum at intermediate scales between the Rossby 
deformation radius [in the ocean typically ~ (10 — 50) km] and the dissipative 
scales (much smaller than 1 km). A thorough discussion about submesoscale 
processes and dynamics can be found in Thomas et al. (2008). Recent high- 
resolution 3D simulations of upper layer turbulence (Capet et al. 2008a, b; Klein 
et al. 2008) have shown that the direct cascade energy spectrum flattens from 
fc~ 3 to k~ 2 for order O(l) Rossby number R a = U/fL, where U is the typical 
velocity difference on a characteristic length L of the flow and / ~ 10~ 4 s _1 is 
the Coriolis parameter. 

Our main purpose is to exploit the MONDO drifter trajectories, shown in 
Fig. [21 to examine relative dispersion by means of several indicators, and discuss 
the consistency of our data analysis in comparison with classical turbulence 
theory predictions, model simulations and previous drifter studies available for 
different regions of the ocean. 

This paper is organized as follows: in section [5] we recall the definitions of 
the major indicators of the Lagrangian dispersion process; in section[3]we give a 
description of the MONDO drifter Lagrangian data; in section 2] the outcome of 
the data analysis is presented; in section [3J we report the analysis of the ocean 
model Lagrangian simulations in comparison with the observative data; and, 
in section El we outline and discuss the main results we have obtained in the 
present work. 

2 Lagrangian dispersion indicators 
2.1 One-particle statistics 

Let r = (x, y, z) be the position vector of a Lagrangian particle, in a 3D space, 
evolving according to the equation dr/dt = U(x, y, z, t), where \J(x,y, z,t) is a 
3D Eulerian velocity field, and let us indicate with v(£) the Lagrangian velocity 
along the trajectory r(t). Let us imagine, then, a large ensemble of N 3> 1 
Lagrangian particles, passively advected by the given velocity field, and refer, 
for every statistically averaged quantity, to the mean over the ensemble. 

The autocorrelation function of a Lagrangian velocity component v can be 
defined, for t\ > t , as: 

C(h,t ) = [(v(h)v(to)) - MM) 2 ] / [M*o) 2 > - (v(t )) 2 ] (1) 

In case of stationary statistics, C(ti, to) depends only on the time lag t = t\ — to- 
The integral Lagrangian time tj, is the time scale after which the autocorrelation 
has nearly relaxed to zero. Typically tl can be estimated as the time of the 
first zero crossing of C(ti, to), or, alternatively, as the time after which \C(t\, t )\ 
remains smaller than a given threshold. 
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Absolute dispersion can be defined as the variance of the particle displace- 
ment relatively to the mean position at time i: 

(A 2 (t)) = {[r(t)-r(0)] 2 )-{[r(t)-r(0)}) 2 (2) 

In the limit of very small times, absolute dispersion is expected to behave as 
follows: 

(A(t) 2 ) ~ alt 2 , for t « t l (3) 

where tl is the Lagrangian autocorrelation time and a\ is the total Lagrangian 
velocity variance. The ballistic regime ([3]) lasts as long as the trajectories save 
some memory of their initial conditions. In the opposite limit of very large 
times, when the autocorrelations have relaxed to zero and the memory of the 
initial conditions is lost, we have: 

(A(t) 2 ) ~ 21ft, for t > t l (4) 

where K is the absolute diffusion coefficient (Taylor 1921). 

Although single particle statistics give information about the advective trans- 
port, mostly because of the largest and most energetic scales of motion, two (or 
more) particle statistics give information about the physical mechanisms acting 
at any scale of motion, compatibly with the available resolution. 

2.2 Two-particle statistics 

Let us indicate with R(t) = \\r^(t) — r^ 2 ^(t)|| the distance between two trajec- 
tories at time t. Relative dispersion is defined as the second order moment of 
R(t): 

( J R 2 (t)) = (||r( 1 )(t)-r( 2 )(t)|| 2 ), (5) 

where the average is over all the available trajectory pairs (r^,!^ 2 )). 

In the small scale range, the velocity field between two sufficiently close tra- 
jectories is reasonably assumed to vary smoothly. This means that, in nonlinear 
flows, the particle pair separation typically evolves following an exponential law: 

(R 2 {t)) ~ e i(2)t (6) 

where, from the theory of dynamical systems, L(2) is the Generalized Lyapunov 
Exponent of order 2 (Bohr et al. 1998). When fluctuations of the finite-time 
exponential growth rate around its mean value are weak, one has L(2) ~ 2Al, 
where Al is the (Lagrangian) Maximum Lyapunov exponent (MLE; Boffetta et 
al. 2000). Notice that for ergodic trajectory evolutions the Lyapunov expo- 
nents do not depend on the initial conditions. If Ai > (except for a set of 
zero probability measure) we speak of Lagrangian chaos. The chaotic regime 
holds as long as the trajectory separation remains sufficiently smaller than the 
characteristic scales of motion. 
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In the opposite limit of large particle separations, when two trajectories are 
sufficiently distant from each other to be considered uncorrelated, the mean 
square relative displacement behaves as: 

(R 2 {t)) ~4K E t, forf-»oo (7) 

where we indicate with Ke the asymptotic eddy-diffusion coefficient (Richard- 
son 1926). 

At any time t, the diffusivity K{t) can be defined as: 

with K(t) -» K E for t -> oo. 

If the velocity field is characterized by several scales of motion, relative 
dispersion in the intermediate range, i.e. between the smallest and the largest 
characteristic length, depends on the properties of local velocity differences, 
i.e. the mean gradients on finite scale. For instance, in 3D fully developed 
turbulence (Frisch 1995), relative dispersion follows the so-called Richardson's 
law: 

(#(*)> ~* 7 , (9) 

with 7 = 3, as long as the trajectory separation lies in the inertial range of 
the energy cascade (Richardson 1926) from large to small scales. It is worth 
to remark that Richardson's law also holds in the inverse cascade range (from 
small to large scales) of 2D turbulence because, in that case as well, the energy 
spectrum follows Kolmogorov's A: -5 / 3 scaling, exactly as in the inertial range 
of 3D turbulence (Kraichnan 1967). Any power law of the type J9| for (R 2 (t)) 
with 7 > 1 is known as super-diffusion. 

2.3 Scale-dependent indicators 

The Finite-Scale Lyapunov Exponent (FSLE) has been formerly introduced as 
the generalization of the MLE A for non-infinitesimal perturbations (Aurell et 
al. 1996). If S is the size of the perturbation on a trajectory in the phase space 
of a system, and (r(<5)) is the phase space averaged time that S takes to be 
amplified by a factor p > 1, then the FSLE is defined as 

The quantity t(S) is the exit time of the perturbation size from the scale S, 
and it is defined as the first arrival time to the scale p ■ S, with p ~ 0(1). The 
computation of the expectation value of the growth rate at a fixed scale, which 
justifies the definition (|10p . is described in Boffetta et al. (2000). As far as 
Lagrangian dynamics are concerned, the evolution equations of the Lagrangian 
trajectories form a dynamical system whose phase space is the physical space 
spanned by the trajectories. In this context, the analysis of relative dispersion 



6 



can be treated as a problem of finite-size perturbation evolution, with scale- 
dependent growth rate measured by the FSLE. The first who had the idea to 
measure the relative dispersion, or, equivalently, the diffusivity, as a function of 
the trajectory separation was Richardson (1926). The FSLE is fundamentally 
based on the same principle. Recently, the use of fixed-time and fixed-scale 
averaged indicators of relative dispersion in various contexts, from dynamical 
systems to observative data in ocean and atmosphere, have been reviewed and 
discussed in several works (Artale et al. 1997; Boffetta et al. 2000; Lacorata et 
al. 2001, 2004; LaCasce and Ohlmann 2003, LaCasce 2008). By a dimensional 
argument, if relative dispersion follows a (R 2 (t)) ~ t 2 ^ scaling law, then the 
FSLE is expected to scale as X(S) ~ . For example, in the case of stan- 
dard diffusion we expect (3 — 2; in Richardson's super-diffusion, j3 — 2/3; in 
ballistic or shear dispersion we have /3 = 1. Chaotic advection means expo- 
nential separation between trajectories. In terms of FSLE this amounts to a 
scale-independent X(S) =constant: that is, (3 — »• 0. In the limit of infinitesimal 
separation, the FSLE is nothing but the MLE, i.e. X(5) ~ Xl (Aurell et al. 
1996). Under these conditions, relative dispersion is said to be a non-local pro- 
cess, since it is determined by velocity field structures with a characteristic scale 
much larger than the particle separation. On the contrary, when the growth of 
the distance between two particles is mainly driven by velocity field structures 
of the same scale as the particle separation, relative dispersion is said to be a 
local process. The super-diffusive processes occurring in 2D and 3D turbulence 
are phenomena of this type. 

An indicator related to the FSLE is the mean square velocity difference 
between two trajectories as function of their separation. Indicating with r^, 
r ( 2 ) ( v ( 1 ) j v ( 2 ) the positions and the Lagrangian velocities, respectively, of two 
particles 1 and 2 at a given time, we define the Finite-Scale Relative Velocity 
(FSRV) at scale 5: 

([Ay(i)f) = ([v (1) -v< 2 )] 2 ) (ii) 

where the average is over all trajectory pairs fulfilling the condition R(t) = 
||rW(£) — r( 2 )(£)|| = 5 at some time t. From the FSRV a scale-dependent 
diffusivity can be formed as K{5) = (1/2) ■ 5 ■ ([AV(5)} 2 ) 1 ^ 2 and compared to 
the classical time-dependent diffusivity K(t) defined in ([5]). 

We would like to end this section with a remark. It is physically reason- 
able to assume that, for statistically homogeneous flows, the behavior of the 
relative dispersion changes in correspondence of certain characteristic lengths 
of the system, rather than after certain time intervals. The major shortcoming 
when measuring the mean square relative displacement at fixed time consists 
in averaging together particle pairs that may have very different separations, 
i.e. may belong to very different dispersion regimes. On the contrary, the mean 
dispersion rate at fixed scale, or other scale-dependent related quantities, gener- 
ally have the benefit of being weakly affected by possible overlap effects. This is 
of particular relevance since the presence of scaling laws in the pre- asymptotic 
regime raises questions about the characteristics (e.g., the kinetic energy spec- 
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trum) of the underlying velocity field. 



3 Data set and local oceanography 

The MONDO Project is a campaign planned by PROOCEANO for ENI Oil do 
Brasil to provide unrestricted access to oceanographic data and subsidize not 
only the companies operational goals but also the scientific community. The 
Lagrangian experiment consisted in the deployment of 39 WOCE SVP drifters 
(Sybrandy and Niiler 1991) in the surroundings of an oil drilling operation (24° 
24' 50" S and 43° 46' 22" W) in proximity of the BC. The drifters are of 
the holey-sock type, composed by a surface floater, a tether and a submerged 
drogue, following the design criteria proposed by Sybrandy and Niiler (1991). 
The drogue length is 6.44 m, centered at 15 m deep, in order to represent the 
first 20 m current information. The drifters are equipped with a GPS device 
and IRIDIUM transmission, providing a positioning accuracy of 7 meters and an 
acquiring rate of 3 hours. From the 39 drifters, 14 were deployed individually 
in a 3-day frequency and 25 were deployed in clusters of 5, every 12 days. 
Justification for the 3-day frequency deployment is based upon the 6.5 days 
energy peak in the wind variability in the region (Stech and Lorenzetti 1992). 
Data used in this study comprehend the period within 21 September 2007 - 
14 November 2008 and passed the quality control proposed by Hansen and 
Poulain (1996) in the Global Drifter Program (GDP) database. High frequency 
components have been removed by applying a Blackman low-pass filter with a 
window of 15 points (45 hours). 

The domain explored by the MONDO Project drifters mainly corresponds 
to that of the BC and to the area where the southward flowing BC meets the 
northward flowing MC, forming the BMC, which results in an eastward current 
feeding the SAC. BC is a warm western boundary current, flowing southward 
and meandering over the 200 m isobath (Peterson and Stramma 1991; Lima et 
al. 1996). As reported by de Castro (2006) and Garfield (1990) the velocities 
inside the BC have values between 25 and 80 cm s _1 . BC flux increases with 
latitude (Assireu 2003; Miiller et al. 1998; Gordon and Greengrove 1986) and 
veers towards east when meeting the colder Malvinas Current around 40°S. 
The BMC is a highly energetic zone with strong thermal gradient and great 
variability, playing an important role in weather and climate of South America 
(Legeckis and Gordon 1982; Pezzi et al. 2005; Piola et al. 2008). Both BC and 
BMC show intense mesoscale activity with eddies detaching from both sides of 
the flow. Lentini et al. (2002) used satellite derived Sea Surface Temperature 
(SST) to estimate an average of 7 rings per year, with lifetimes ranging from 
11 to 95 days with major and minor radii of 126 ± 50 km and 65 ± 22 km 
respectively. Assireu (2003) estimates a diffusion coefficient between 6 • 10 6 
and 9.1 • 10 7 cm 2 s _1 , a Lagrangian time scale between 1 and 5 days and a 
characteristic (Eulerian) length scale varying from 19 to 42 km, which agrees 
with the meridional variation of the 1 st internal Rossby radius of deformation 
in the region (Houry et al. 1987). Studies addressing the energetics in BC and 
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Figure 3: Mass center coordinates in longitude (9) and latitude ((f)), averaged 
over all drifters, as function of the time lag since the release. The time sampling 
is At —1/8 day. Error bars are the standard deviations on the mean values. 



BMC flows report great spatial and temporal variability (see, e.g., Assireu et al. 
2003; Piola et al. 1987; Stevenson 1996) with EKE varying from 200 to 1200 
cm 2 s~ 2 and representing 70-90% of the Total Kinetic Energy (TKE). Recently, 
Oliveira et al. (2009) extended the analysis made by Piola et al. (1987) and 
suggested that EKE is usually lower than Mean Kinetic Energy (MKE) in the 
BC and always greater in the BMC. 



4 Data analysis 

4.1 Single particle statistics 

From a pre-analysis we eliminated two drifters that very soon lost their drogue, 
so the set of active drifters is reduced to 37 units. The Lagrangian velocity statis- 
tics, reconstructed along the drifter trajectories by means of a finite difference 
scheme, provides standard deviations, in the zonal and meridional directions, of 
order ~ 25 cm s , compatible with estimates obtained from other observative 
data available for the same ocean region (Figueroa and Olson 1989; Schafer and 
Krauss 1995). 

The drifters initially tend to move southward, driven by the BC; then they 
are nearly stopped by the BMC, and, last, they tend to be transported east- 
ward by the SAC. The time evolution of the drifter mean position coordinates 
in longitude and latitude is shown in Fig. [3] During the first phase of dispersion, 
the mean meridional displacement grows almost linearly southward (BC), while 
the mean zonal displacement is almost constant; after about 120 days from the 
release time, the mean meridional displacement tends to saturate (BMC), while 
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Figure 4: Lagrangian velocity autocorrelation function for zonal and meridional 
components. The time sampling is At = 1/8 day. 

the mean zonal displacement tends to grow almost linearly eastward (SAC). 
Given the differences between BC, BMC and SAC in terms of energy and cir- 
culation patterns, and, in general, between the boundary of the Subtropical 
Gyre and its interior, the drifters do not sample a statistically homogeneous 
and stationary domain. In order to evaluate the Lagrangian autocorrelations, 
we divide the trajectory lifetime into time windows of 10 days width. Then we 
compute ([TJ) where mean and variance of the Lagrangian velocity components 
are recalculated for each window, and then we take the average, for a fixed time 
lag t = t\ — in, over all the windows for all trajectories. The resulting functions 
are plotted in Fig. S) We notice that, since most of the statistics regards trajec- 
tory segments within or near the BC, the motion is slightly more autocorrelated 
in the meridional than in the zonal direction. However, the order of the integral 
time scale is the same for both components. Without discussing the various 
estimates that can be evaluated using different criteria, we identify the two in- 
tegral time scales around a value tl ~ 5 days, with a difference between zonal 
and meridional components of order 0(1) days. This implies that, for time lags 
significantly larger than tl , the Lagrangian velocities along a trajectory become 
practically uncorrelated. 

Absolute dispersion, i.e. the mean square fluctuation around the mean po- 
sition, for zonal and meridional coordinates is reported in Fig. [5] 

In the present work we are not interested in reproducing estimates of small 
scale diffusivity, obtained by reconstructing a mean velocity field, in general not 
spatially uniform, and calculating the turbulent components from the difference 
between the drifter local velocities and the local mean field. Wherever we speak 
of diffusion coefficients, referring to ([2] ) or ([5] ) . we mean the effective diffusion 
coefficients of order ~ 0(UL) where L and U are determined, respectively, 
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Figure 5: Mean square drifter displacement from the mass center position, 
(A 2 (t)), in zonal and meridional components. The time sampling is At — 
1/8 day. 



by the characteristic size and rotation velocity of the largest eddies. Absolute 
dispersion, being dominated by the large scale features of the velocity field, 
reflects the anisotropy of the currents encountered by the drifters. The ballistic 
dispersion ~ t 2 scaling is plotted as reference for autocorrelated motions on time 
lags smaller than 10 days. The diffusion ~ t scaling is a good approximation for 
meridional dispersion on time lags within 10 — 60 days, but it is not good for 
zonal dispersion on any time lag. After about 60 days, meridional dispersion 
tends to saturate, while zonal dispersion starts growing again as ~ i 2 , likely 
because, on those time lags, on average, the drifters have reached the BMC, 
have stopped their southward meridional run along the BC and start moving 
eastward along the SAC, the natural consequence of which is the saturation of 
the meridional component and a second ballistic (or shear) regime for the zonal 
component. 

4.2 Two-particle statistics 

Relative dispersion for three different initial separations is shown in Fig. [51 
Distances between two points on the ocean surface are calculated as great circle 
arcs, according to the spherical geometry standard approximation. A way to 
increase relative dispersion statistics is to measure flSJ) also for so-called chance 
pairs, i.e., particle pairs which may be initially very distant from each other but 
that, following the flow, happen to get sufficiently close to each other at some 
later time after their release, in some random point of the domain (LaCasce 
2008). The maximum number of pairs counted in the statistics depends on the 
initial threshold: 24 pairs for R(0) < 1 km, 30 pairs for R(0) < 5 km and 39 
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Figure 6: Drifter relative dispersion (R 2 (t)) for initial separations i?(0) < 1 km, 
i?(0) < 5 km and i?(0) < 10 km. The quantity L(2) ~ 1.2 day -1 corresponds 
to \l ~ 0.6 day -1 . The time sampling is At — 1/8 day. Error bars are the 
standard deviations on the mean values. 

pairs for R(0) < 10 km. The dependence of (R 2 (t)) on R(0) is well evident. 
The early regime of relative dispersion is reported in Fig. [7j The value of the 
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Figure 7: Early regime of drifter relative dispersion (R 2 (t)), for initial separa- 
tions i?(0) < 1, < 5 and i?(0) < 10 km, with L(2) ~ 2A L ^ 1.2 day" 1 . The 
time sampling is At = 1/8 day. Error bars are the standard deviations on the 
mean values. 
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Figure 8: Drifter relative diffusivity K{t) in zonal and meridional components, 
for initial separations R(0) < 5 km. The t 2 scaling is plotted as reference to the 
Richardson regime. The time sampling is At = 1/8 day. 



MLE Al — 0.6 day -1 is of the same order as other estimates obtained in various 
ocean regions, as discussed later. 

The exponential growth has been adapted to the R(0) < 5 km curve which 
has a fair compromise between having a sufficiently small initial separation and, 
at the same time, an acceptable pair statistics. 

Relative diffusivity in the zonal and meridional directions, for initial sepa- 
rations R(0) < 5 km, is plotted in Fig. [5J Although the curves do not display 
a clean scaling behavior, in the time range 10 — 100 days the diffusivities show 
some resemblance with the t 2 law expected in the Richardson dispersion regime. 
We will examine now the results obtained with scale-dependent indicators. 

The FSLE has been evaluated for the same initial thresholds 1, 5 and 10 km; 
the results are reported in Fig. [S] Particle pairs are counted in the statistics 
if, at a given time, they happen to to be separated by a distance smaller than 
the fixed threshold. The density of scales r5 is fixed by setting p = \[2. The 
maximum number of pairs considered varies with the initial threshold: 24 pairs 
for R(0) < 1 km, 30 pairs for R(0) < 5 km and 39 pairs for R(0) < 10 km. 
Larger initial thresholds (viz., > 100 km) could increase the statistics of an order 
~ O(10), but we will not consider this opportunity since, with too large initial 
thresholds, mesoscale dispersion rates may be affected by the inhomogeneity of 
the flow. In other words, at a certain separation scale two drifters might be still 
spatially correlated in some regions of the domain, but totally uncorrelated in 
some other ones, depending on the structure of the flow. The parameter p cannot 
be chosen arbitrarily close to 1 for obvious limitations because of finite resolution 
of the data. We have checked that, with p = y/2, "clipping events", i.e., suspect 
transitions between two neighbouring scales in only one time step At = 1/8 days, 
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Figure 9: FSLE for i?(0) < 1, < 5 and < 10 km and amplification factor 
p = \/~2. The small scale plateau corresponds to exponential separation with 
rate \l — 0.6 day -1 . The <~ <5~ 2 / 3 scaling corresponds to the Richardson law. 
The S" 1 ^ 2 scaling corresponds to a k~ 2 spectrum. Error bars are the standard 
deviations on the mean values. 

are not significant. Taking a value 5r — 30 km for the Rossby radius, we can see 
that mesoscale dispersion (scales > 5r) is compatible with a S~ 2 ^ 3 Richardson 
regime, as it occurs in the inverse cascade range of QG turbulence. Cut-off due 
to finite lifetime of the trajectories occurs before reaching <~ 10 3 km separations. 
The range (1 — 100) km is characterized by a step-like behavior of the FSLE. The 
first plateau in the submesoscale range ~ 0(1) km corresponds to exponential 
separation with A ~ 0.6 day -1 , in agreement with the result obtained before for 
(R 2 (t)). The second plateau at ~ 0.15 day -1 on ~ O(Sr) scales is compatible 
with exponential separation due to a k~ 3 direct cascade of QG turbulence. It 
is unclear why there exist two quasi-constant levels for A(<5), which produce a 
significant variation of the exponential growth rate when S decreases from meso 
to submesoscales. Possible interpretations of this behavior may be found in 
the action of submesoscale vortices, even though the FSLE does not display 
a continuous cascade connecting these two scale ranges. The ~ 5~ 1 ' 2 scaling 
is plotted as a reference to a k~ 2 energy spectrum. Another remark concerns 
the gap at 6 <~ 100 km between the end of the plateau and the beginning of 
the ~ <5~ 2 / 3 regime. Inverse cascade implies that smaller eddies, formed by 
instability at the Rossby scale, tend to merge together in larger eddies up to 
scales in the range (10 2 — 10 3 ) km. Assuming this actually occurs, we also have 
to consider that the BC and BMC generate rings of size ~ O(10 2 ) km which 
detach from the main currents. Such mesoscale rings can have trapping effects 
on the drifters, at least to some extent, which would determine an abrupt fall 
of the dispersion rate (FSLE). The presence of the BMC can also have a role 
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Figure 10: FSRV computed together with the FSLE. The FSRV is the mean exit 
velocity of the distance between two drifters from the "shell" 5. The scalings 
S 2 , <5 2 / 3 and S correspond to k~ 3 , fc~ 5 / 3 and k~ 2 energy spectra, respectively. 
Error bars are the standard deviations on the mean values. 



in temporarily inhibiting the growth of the two particle separation, until the 
drifters find their way eastward along the SAC. 

Let us now consider the finite-scale (mean square) velocity difference, Fig. 
1101 The Lagrangian velocity components are reconstructed from the trajectories 
by means of a finite difference scheme. We notice, first, that the FSRV returns 
back a slightly cleaner picture than other indicators. This may depend, for 
instance, on the fact that the amplification time of the separation between 
two drifters passing from S to p ■ 5 is affected by large fluctuations due to the 
non homogeneous characteristics of the flow, while the (mean square) relative 
velocity at scale S is less affected by the history of the dispersion between two 
scales. The two "valleys" at nearly 10 and 100 km leave room to hypothesis 
of trapping events by structures of size comparable with those scales. The 
FSRV is consistent with Richardson dispersion on scales larger than the Rossby 
radius (except for the "valley" at 100 km); with exponential separation for 
two subranges ~ (10 — 50) km and ~ (1 — 5) km; with a 5 scaling in the 
subrange ~ (5 — 10) km, in turn consistent with a k~ 2 energy spectrum. The 
~ (10 — 50) km exponential regime can be associated to a fc~ 3 direct cascade, 
while the ~ (1 — 5) km exponential regime can be associated to the action of 
submesoscale vortices. From the FSRV an "equivalent Lagrangian spectrum" 
Ez,{k) = ([AV(k)] 2 ) /k can be defined, by dimensional arguments, replacing S 
with 2w/k, Fig. [TT] The same scenario formerly indicated by the FSRV is 
reproduced, in k space, by E^k) as well. 

We compare now the diffusivity (see Fig. [T2|) computed in both ways: as a 
fixed-time average (|8]) and as a fixed-scale average from the FSRV. Both quan- 
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Figure 11: Equivalent Lagrangian spectrum denned from the FSRV. The inverse 
cascade (fc~ 5 / 3 ), direct cascade (k~ 3 ) and the submesoscale (fc -2 ) spectra are 
plotted as reference. The Rossby radius ~ 30 km corresponds to a wavenumber 
k ~ 0.2. 
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Figure 12: Diffusivity as a function of the separation: fixed-time average K(t) 
vs 8 = (i? 2 (i)) 1//2 and fixed-scale average K(5) vs S. The <5 4 / 3 and i5 2 correspond 
to fc -5 / 3 and /c -3 spectra, respectively. The <5 3 / 2 scaling corresponds to a k~ 2 
spectrum. 



tities are plotted as functions of the separation between two drifters: K(5) = 
(1/2) • 8 ■ ([AV(S)] ) 1 / 2 with S as the independent variable, and K(t) versus 
S = {R 2 ^)) 1 / 2 where the independent variable is the time t. The diffusivity is 
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less affected by distortions present in the signal than other indicators. The two 
regimes ~ <5 4 / 3 at mesoscales and ~ 5 2 at scales ~ O(Sn) are consistent with the 
inverse cascade fc~ 5 / 3 and direct cascade fc -3 , respectively, as predicted by QG 
turbulence theory. A scaling ~ S 3 / 2 , corresponding to a fc -2 spectrum, connects 
scales of the order of the Rossby radius to the submesoscale, and a scaling ~ S 2 
is the signature of exponential separation driven by velocity field features on 
scales ~ 0(1) km. 

5 Numerical simulations 

To check the consistency of the results obtained from the data analyis of the 
MONDO drifters, we have performed two numerical experiments in which ~ 
O(10 2 ) virtual drifters are "deployed" in a surface current field generated by a 
global operational forecast system using the Hybrid Coordinates Ocean Model 
(HYCOM; Bleck and Benjamin 1993) and the Navy Coupled Ocean Data As- 
similation (NCODA; Cummings 2005). The model grid step is 1/12° (approx- 
imately 7 km) and outputs are available with a 1 day time step. From the 32 
layer global grid, the surface layer is extracted in the area within 45 o -20°S' and 
60°-30°W for the same period of drifters trajectories (20 September 2007 - 21 
October 2008). Conversion from Earth coordinates to meters is accomplished 
following the great circle distance approximation. The integration algorithm of 
the Lagrangian trajectories uses a fixed timestep At = 1/24 day and tri- linear 
interpolations. The maximum integration time is 400 days but particles reaching 
the shoreline are stopped and removed from the subsequent integration. 

We indicate the numerical experiments with El and E2. In the first one 
(El) the drifters are uniformly deployed in an area of about (10 x 10) km 2 cen- 
tered around a position corresponding to the mean initial location of MONDO 
drifters. The average initial distance between particle pairs is (R(0)) — 5 km. 
The lifetime of trajectories is between 150 days and 200 days. In the second 
experiment (E2) the initial distribution of the drifters is characterized by larger 
separations, namely comparable to the spacing of the numerical grid (~ 10 km); 
the average initial distance between particle pairs is (R(0)) — 40 km and the 
duration of trajectories is 250 — 400 days. 

Here we present some comparisons between the analysis of data from model 
trajectories and those from the actual drifters, focusing on the relative dispersion 
process. First we consider a fixed-time indicator, namely relative dispersion 
(R 2 (t)), for trajectories selected with small enough initial separation, that is for 
i?(0) < 10 km. As it can be seen in Fig. [T3J the model outcome is in reasonable 
agreement with the behaviour of actual drifters. In particular, experiment El 
is capable of reproducing the early growth of relative dispersion but it cannot 
catch the late time behaviour due to too short trajectories. On the other hand, 
with experiment E2 we obtain a nicer agreement at late times, thanks to the 
longer lifetime of virtual trajectories in this case, but we loose the agreement 
at early times because now the initial separation between particle pairs is on 
average much larger. 
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Figure 13: Relative dispersion (R 2 (t)) for initial separations R(0) < 10 km 
for MONDO drifters and virtual drifters from numerical experiments El and 
E2. For virtual drifters errors are of the order of point size; time sampling is 
At = 1/8 day for MONDO drifters and At = 1/24 day for virtual drifters. 

For what concerns fixed-scale indicators, we consider the FSLE and the 
FSRV (Figures [Ml [T5|) . To increase the statistics we now select trajectories 
with a larger initial separation, that is for i?(0) < 50 km (similar results are 
found for smaller values of i?(0), though they are more noisy). The behaviours 
of both FSLE and FSRV support a double cascade scenario on the same scales 
as those found with MONDO drifters. The plateau value of FSLE at scales 
O(Sn) is very close to the one found in the real experiment (X(S) ~ 0.15 day -1 ). 
At larger scales, for both numerical experiments El, E2 the behaviour of FSLE 
is compatible with a scaling law A(<5) ~ <V 2 / 3 supporting an inverse energy 
cascade process. Experiment E2, which is characterized by longer trajectories, 
shows a clearer scaling, thanks to a larger number of pairs reaching this range of 
large scales. Mean square velocity differences in the same range of scales display 
a reasonably clear r5 2 / 3 scaling, also supporting an inverse energy cascade, with 
values close to those found with MONDO drifters. At separations smaller than 
the Rossby deformation radius, both indicators point to the presence of a direct 
enstrophy cascade: the FSLE is constant and the FSRV behaves as ([AV(<5)] 2 ) ~ 
S 2 . This only partially agrees with the results found for real drifters, namely only 
in the scale range 10 km < S < 50 km. At subgrid scales, velocity field features 
are not resolved and relative dispersion is necessarily a nonlocal exponential 
process driven by structures of size of the order of (at least) the Rossby radius. 
Correspondingly, the FSLE computed on model trajectories does not display 
the higher plateau level at scales smaller than 10 km. 

Finally, the behaviour of the relative diffusivity K{8) as a function of the 
separation 8 is shown in Fig. [TC] for both MONDO and virtual drifters. Here 
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Figure 14: FSLE for R(0) < 50 km and amplification factor p = ^2 for MONDO 
drifters and virtual drifters from numerical experiments El and E2. For virtual 
drifters errors are of the order of point size. The large-scale saturation (El) 
depends on the value of the trajectory integration time. 



K(S) is computed from the mean square velocity difference, as described in 
Sections [2}: and [4}d. Model and experimental data again are in agreement at 
scales larger than the numerical space resolution (~ 10 km). Indeed, in both 




Figure 15: FSRV for R(0) < 50 km for MONDO drifters and virtual drifters 
from numerical experiments El and E2. For virtual drifters errors are of the 
order of point size. 
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numerical experiments El and E2 we find scaling behaviours compatible with a 
QG double cascade: K(5) <~ 5 4 / 3 (corresponding to Richardson's super-diffusion 
in an inverse cascade regime) for S > Sr, and K(5) ~ S 2 (corresponding to a 
direct cascade smooth flow) for 6 < 6r. At variance with the outcome of the real 
experiment with MONDO drifters, here we are unable to detect any significant 
deviations from the QG turbulence scenario at small scales. 




Figure 16: Diffusivity K{8) as function of the separation 5 for MONDO drifters 
and virtual drifters from numerical experiments El and E2; here R(0) < 50 km. 

The failure of the model to reproduce flow features at very small scale is 
of course due to its finite spatial resolution, which is of the order of 10 km. 
Below this length scale, the velocity field computed in the model is smooth, 
while the one measured in the real ocean clearly displays active scales also in 
the range (1 — 10) km. Nevertheless, the overall conclusion we can draw from the 
above comparisons is that the characteristics of the relative dispersion process 
found with MONDO drifters agree with those obtained with an Ocean General 
Circulation Model (OGCM) for scales <5 > 10 km. 



6 Discussion and conclusions 

Lagrangian dispersion properties of drifters launched in the southwestern cor- 
ner of the South Atlantic Subtropical Gyre have been analyzed through the 
computation of time-dependent and scale-dependent indicators. 

The data come from a set of 37 WOCE SVP drifters deployed in the Brazil 
Current around (24°S, 44°W) during Monitoring by Ocean Drifters (MONDO) 
Project, an oceanographic campaign planned by PROOCEANO and supported 
by ENI Oil do Brasil in relationship with an oil drilling operation. The ex- 
perimental strategy of deploying part of the drifters in 5-element clusters, with 
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initial separation between drifters smaller than 1 km, allows to study relative 
dispersion on a wide range of scales. 

Single-particle analysis has been performed by computing classic quantities 
like Lagrangian autocorrelation functions and absolute dispersion, defined as 
the variance around the drifter mean position, as a function of the time lag 
from the release. Velocity variances (<~ 500 cm 2 s~ 2 ) and integral time scales 
(tl ~ 5 days) are compatible with the estimates obtained in the analysis of the 
FGGE drifters (Figueroa and Olson 1989; Schafer and Krauss 1995). Anisotropy 
of the flow is reflected in the different behavior of the zonal and meridional 
components of the absolute dispersion. Being the MONDO drifters advected 
mostly by the boundary currents surrounding the Subtropical Gyre, the Brazil 
Current first and the South Atlantic Current later, the meridional component of 
the absolute dispersion is dominant as long as the mean drifter direction is nearly 
southward (BC), while the zonal component dominates (see the appearance of 
the late ballistic regime) as the mean drifter direction is nearly eastward (SAC). 
Early time advection is modulated, also, by the response of the currents to a 
wind forcing of period <~ 6.5 days, a characteristic meteorological feature of the 
BC dynamics. 

Two-particle analysis has been performed by means of both fixed-time and 
fixed-scale averaged quantities. Classic indicators like the mean square rela- 
tive displacement and the relative diffusivities between two drifters as functions 
of the time lag from the release give loose information about early phase ex- 
ponential separation, characterized by a mean rate Al — 0.6 day -1 , and long 
time dispersion approximated, to some extent, by Richardson super-diffusion 
before the cut-off due to the finite lifetime of the trajectories. Evidence of a 
small scale exponential regime for relative dispersion is common to other drifter 
studies for different ocean regions (LaCasce and Ohlmann 2003; Ollitrault et 
al. 2005, Koszalka et al. 2009, Lumpkin and Elipot 2010). Scale-dependent 
indicators return back a cleaner picture, compatibly with the limited statistics 
allowed by the experimental data and the non homogeneous and non stationary 
characteristics of the flow. 

The FSLE displays a mesoscale [(100 — 500) km] regime compatible with 
Richardson super-diffusion, Lagrangian counterpart of the 2D inverse cascade 
scenario characterized by a fc~ 5 / 3 energy spectrum. At scales smaller than 100 
km the FSLE has a step-like shape, with a first plateau at level <~ 0.6 day -1 in 
the submesoscale range <~ 0(1) km, and a second plateau at level ~ 0.15 day -1 
for scales comparable with the Rossby radius of deformation (Sr ~ 30 km). 
Constant FSLE in a range of scales corresponds to exponential separation. The 
~ 0(6r) plateau could be related to the 2D direct cascade characterized by a 
fc~ 3 energy spectrum, while the origin of the ~ 0(1) km plateau is likely related 
to the existence of submesoscale features of the velocity field, the role of which 
has been recently assessed by means of high-resolution 3D simulations of upper 
ocean turbulence at Rossby numbers R <~ 0(1) (Capet et al. 2008a, b; Klein 
et al. 2008). The FSLE does not display a clean continuous cascade scaling 
<~ (5 -1 / 2 , corresponding to a k~ 2 energy spectrum, from sub to mesoscales; 
however, it highlights the existence of scales of motion hardly reconcilable with 
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the QG turbulence scenario, as analogously assessed also by Lumpkin and Elipot 
(2010) for drifter dispersion in the North Atlantic. 

The FSRV measures the mean square velocity difference at scale S. The 
scaling of the FSRV is related to the turbulent characteristics of relative dis- 
persion. The FSRV behavior is, under a certain point of view, cleaner than 
that of the FSLE, but it is affected by the presence of two "valleys", roughly 
at ~ 10 km and at ~ 100 km, likely associated to trapping events (the same 
features are present also in the FSLE behavior). Coherent structures on scales 
<~ 10 and <~ 100 km may be responsible of the "fall" in the relative dispersion 
rate. These scales are of the order of the Rossby radius and of the mesoscale 
rings that detach from the current systems, respectively. We must also consider 
the role of the Brazil-Malvinas Confluence which tends to inhibit the growth of 
the dispersion as the drifters, initially flowing southwestward along the Brazil 
Current, encounter the northeastward flowing Malvinas Current before being, 
eventually, transported eastward along the South Atlantic Current. The FSRV 
displays a mesoscale ~ 5 2 ' 3 scaling (except for the "valley"), compatible with a 
£,-5/3 mverse cascade; a ~ S 2 scaling in the two subranges ~ (10 — 50) km and 
<~ (1 — 5) km, which correspond to the step-like shape of the FSLE; a ~ S scaling 
which, to some extent, is compatible with a k~ 2 regime connecting scales of the 
order of the Rossby radius to the submesoscale (below which the velocity field 
becomes approximately smooth). The equivalent Lagrangian spectrum E^k), 
formed by dividing the FSRV (essentially the relative kinetic energy at scale 5) 
by the wavenumber fc = 2-7r/<5, returns back the same scenario in k space. 

Last, we have compared the scale-dependent relative diffusivity K{8) con- 
structed from the FSRV, for which 5 is the independent variable, with the classic 
relative diffusivity K(t) seen as a function of the mean separation between two 
drifters (i? 2 (t)) 1 / 2 , for which t is the independent variable. What emerges from 
the analysis of the diffusivities is that, in the mesoscale range, loosely from 
the Rossby radius up to scales <~ O(10 2 ) km, K(5) ~ <5 4 / 3 , corresponding to 
the /c~ 5 / 3 inverse cascade; in the submesoscale range ~ (1 — 5) km, and in a 
limited subrange from about the Rossby radius down to <~ 10 km, K(S) <~ 5 2 , 
corresponding to exponential separation; in the subrange ~ (5 — 10) km, the 
scaling K{5) <~ S 3 ^ 2 is compatible with a k~ 2 regime. The k~ 2 spectrum is 
a characteristic of upper ocean turbulence when the Rossby number is order 
<~ 0(1), as recently assessed with high-resolution 3D model simulations (Capet 
et al. 2008a, b; Klein et al. 2008), which connects the mesoscale to the sub- 
mesoscale ~ 0(1 — 10) km (McWilliams 1985). Below the submesoscale the 
velocity field is reasonably assumed to vary smoothly. A rough estimate of the 
Rossby number associated to the MONDO drifter dynamics gives a value (at 
least) R Q = U/(Lf) ~ C^IO" 1 ), taking U ~ 0.3 m/s, L ~ 30 km and / - 10~ 4 
1/s, which is, nonetheless, considerably larger than the typical Rossby number 
in open ocean. Although the drifter data analysis does not show a clear evidence 
of a relative dispersion regime corresponding to the k~ 2 spectrum, the presence 
of velocity field features of size comparable to submesoscale vortices is reflected, 
to some extent, by the behavior of the small scale relative dispersion process. 

Numerical simulations of the Lagrangian dynamics have been performed 
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with an OGCM of the South Atlantic (Huntley ct al. 2010). The results con- 
cerning the relative dispersion essentially agree with the data analysis of the 
MONDO drifters, within the limits of the available numerical resolution. In 
particular, two-particle statistical indicators such as the FSLE, the FSRV and 
the scale-dependent relative diffusivity K(S), computed on the trajectories of 
virtual drifters, display the same behaviour found for MONDO drifters for scales 
larger than approximatley 10 km, that is larger than the numerical grid spacing. 
Below this length scale, evidently, the model flow field is smooth, hence relevant 
departures from QG turbulence and the role of submesoscale structures cannot 
be assessed. Further investigations on the modeling of submesoscale processes 
would provide extremely useful in order to make a clearer picture of the small 
scale dynamics of the surface ocean circulation in the region. 
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